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We derive solutions for the twin wall linking two tetragonal variants of the cubic-tetragonal fer- 
roelastic transformation, including for the first time the dilatational and shear energies and strains. 
Our solutions satisfy the compatibility relations exactly and are obtained at all temperatures. They 
require four non-vanishing strains except at the Barsch-Krumhansl temperature Tbk (where only 
the two deviatoric strains are needed). Between the critical temperature and Tbk, material in the 
wall region is dilated, while below Tbk it is compressed. In agreement with experiment and more 
general theory, the twin wall lies in a cubic 110- type plane. We obtain the wall energy numerically 
as a function of temperature and we derive a simple estimate which agrees well with these values. 

PACS numbers: 81.30.Kf, 68.35.-p, 62.20.Dc, 61.70.Ng 



Ferroelastic transformations are diffusionless, fipfitr 
order, shape-changing phase changes in the solid statad'0. 
In cubic-tetragonal (C-T) systems like NbsSn, VaSi, In- 
Tl alloys, Fe-Pd alloys and Ni2MnGa, the cubic unit cell 
elongates (or contracts) along one of three axes to form a 
tetragonal unit cell; below the transition temperature Tc 
there are three possible homogeneous products (variants) 
differing only in oricntaticMi. 

Barsch and KrumhanslB (BK) obtained an analytical 
solution for the twin wall linking two tetragonal variants 
of C-T ferroelastics. The dilatational and shear strains 
vanish identically, and the two remaining strains are func- 
tions of a single coordinate. The wall lies in a cubic 110- 
type plane, in agreement with experiment and as already 
known on more general groundsQ. 

The BK solution is valid, however, at only a single 
temperature (T — Tbk)- At any other T, the dilata- 
tional and shear strains are not zero, and the only known 
method to find the C-T twin wall structure requires solv- 
ing the full three-dimensionjaL (3D) partial differential 
equations (previous attemptsErO at a ID solution omitted 
the non-deviatoric strains). Some numerical solutionsLTLl 
of these equations have been obtained, but no results 
have been given for the wall structure. 

That is, 16 years after BK|3, the wall structure is still 
unknown except at Tbk- 

The following solves this long-standing problem, which 
is of considerable physical interest given the large strains 
and large magnetic- field effectstj in Ni2MnGa. Specifi- 
cally, we present a ID solution for the C-T twin wall at 
all T. Our solutions, which include dilatational and shear 
energies and strains, satisfy the compatibility relations by 
virtue of analytical relations that we derive between the 
strains; these relations allow us to reduce the problem to 
the solution of ordinary, rather than partial, differential 
equations. We recover the BK solution at Tbk and we 
present results at both higher and lower T. 

The three paragraphs immediately following define the 
strains and the two parts of the free-energy density. The 
next two paragraphs obtain the key new results of our 
analysis, namely two relations between the strains, and 



the resulting Euler-Lagrange equations. The remain- 
ing paragraphs discuss the results from solution of these 
equations. We find that the twin-wall region is dilated 
near T^ and compressed below Tbk- 

Six strains are required to describe C-T ferroelastics, 
the dilatational strain ei, the deviatoric strains 62 and 
63, and the shear strains 64, 65 and eg. With the coordi- 
nate axes along the four-fold axes, and in the small-strain 
approximation, these are 

ei = (Ul,l + W2,2 + "3,3) /\/3 , 62 = (Uij - ^2,2) / \/2 , 
63 = + U2,2 - 21*3,3) /V6 , eg = (Ui,2 + U2,l) /2 

(1) 

plus obvious expressions for 64 and 65. Here u = 
(ui, 1*2, U3) is the displacement of the material point orig- 
inally at X, and Uij — djUi = dui/dxj. We need also the 
components UJ3 = (wi.2 — U2,i)/2, etc. of the local rota- 
tion u;. 

The free energy F is the integral F = Jy T (fix oi 
the free-energy density over the undeformed volume 
V. For proper ferroelastics (where the strain is the pri- 
mary order parameter), is the sum of strain and strain- 
gradient parts. The strain part is 




In Voigt notation, the dilatational and shear constants 
are Ai = Cn -I- C12 and A4 = 4C44 respectively, with 
both > for stability; the corresponding APWis in the 
density were omitted in previous treatmentSLfuij. The co- 
efficient A2 depends on temperature as A2 — A2{T — Tq), 
where To is the stability limit of the cubic phase and A'2 is 
a material-dependent constant; above Tc, A2 = Cn— C12. 
For A2 > JB2/C2, the energy has only the cubic mini- 
mum (all strains zero). For A2 < JB2/C2, there are in 
addition three degenerate minima symmetrically located 
in the plane of the deviatoric strains: 
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62 = , 63 = 630 
62 = -\/363o/2 , 63 = -630/2 
62 = \/363o/2 , 63 = -630/2 



with 61 = 64 = 65 = 6g 
630 = \B2 + (B^, 



and 
- 4A2C2) 



/(2C2 



(3a) 
(3b) 
(3c) 



(4) 



the tetragonal four-fold axes are in the 3, 1 and 2 di- 
rections respectively. The C-T transition, which is first- 
order, occurs at A2 = |S|/C2 where 630 — |i32/C2- 

The free-energy density requires also strain-gradient 
terms, so that energy is required to introduce variant- 
variant walls (otherwise the system can subdivide into 
arbitrarily fine variants). We keep only the two invariants 
quadratic in the deviatoric strain derivatives: 



(eb)' + K2)' + (e2.3)' 



«i)' + «2)' + (e3,3)' 



(5) 



where 62, 62, 63 and 63 are obtained from 62 and 63 by 
27r/3 rotations about the cubic 111 axis: 



62 = ("2,2 - "3,3) /\/2 = (-62 + \/363) 

62 = ("3,3 - Wl,l) /\/2 = (-62 - \/363) 

("2,2 + U3,3 - 2ui,i) /\/6 = (-63 - \/362) 

("3,3 + "1,1 - 2^2,2) /\/6 = (-63 + \/362) 



/2 
/2 
/2 
/2 



(6) 



Both terms are transparently invariant and non-negative, 
and so we have the stability requirements d2 > and 
^3 > (which diffeii.frQm those in Ref. contact with 
previous treatmentsBflij is made by writing d2 ~ .92 — 93 
and (i3 = 32 + 53 , resulting in 



^.<, = |g2[(V62)' + (V63)^ 



+ ig3{i[(e2,i)'-(63,i)^ 



+ |[(62,2)'-(63,2)']-[(62,3)'-(e3,3)'] 

+ \/3(62, 163a -62,263,2)} • (7) 

We seek the solution linking the variants with four-fold 
axes in the 1 and 2 directions, Eqs. (3b) and (3c); the 
results for this pair are simpler than for the others. The 
method is easily extended to treat a twin band. We as- 
sume a solution with 64, 65, lui and UJ2 — all identically 
zero, and with the other four strains (and lu^) indepen- 
dent of X3. The strains are not independent for phys- 
ical settings like a twin wall where the strains depend 
on position; rather, they are linked by the compatibil- 
ity relations (necessary and sufficient conditions that the 
strains be derivable from the displacement). The nine 
first-order relations, of the form Uijk — Ui^kj, involve 



the first derivatives of the strains and the rotation com- 
ponents uji. The more familiar second-order relations, 
which involve the second derivatives of the strains, are 
easily obtained by differentiation to eliminate the uji . By 
virtue of the above assumptions, only the relations 



d2 (\/26i 

di (V261 



\/362 -f- 

V3e2 + 
di {ei 
d2 (ei 



63) /V6 
63) /V6 

-\/263) 
-\/263 



= dl (66 
= d2 (66 
= , 
= , 



- W3) 



(8) 



need be considered. We try for functions of X = 
xi cos P + X2 sin (3 alone; a similar ID solution is not pos- 
sible for a cubic- tetragonal soliton. The compatibility 
relations are then 

sin (3 (y2ei + V3e2 + 63^ /\/6 - cos /3 (66 + ^3) = Ki , 
cos (3 ^ V261 - V362 + €3^ /VG- sin P {ee - lu^) ^ K2 , 

61 - V263 = K3 . (9) 

The constants Ki, K2 and are evaluated from 
the boundary conditions at X = ±00, namely 62 — 
±iV363o, 63 = -5630, 61 = 66 = 0, and uo^ = ±Vt. One 
finds easily that a solution is possible only if cos 2/3 = 0, 
and so X = (xi±a:2)/\/2; that is, the walls lie in the 110 
or 110 planes. In this way, we find the key new results 
relating 61 and 66 to the deviatoric strains: 



6i(X) = \/2[63(X) + 630/2] , 
66(X) = ±V372[63(X)+ 630/2] , 



(10) 



plus ijJ3{X) = ±e2{X)/\/2. These results are indepen- 
dent of the details of the free-energy density (they apply 
whether or not the dilatational and shear energies appear 
in T). 

Since the compatibility relations are satisfied, we can 
use the density 



^ _ Ads / e^oY A2 
^t™„- — (63 + — j + — 
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(11) 



and minimize F with respect to 62 and 63. Here 



2A1 + IA4, D2 = 1^2 + 1^3, and D3 



jds; adding 



the diagonal invariants ^di{S/ei)^ and id4[(V64)^ 4- 

(V65)^ -I- (V66)^] (with dl and ^4 > 0) to the density 
of Eq. (5) leaves D2 unchanged while adding 2di + |d4 
to D3. The corresponding Euler-Lagrange equations are 



^262 -I- 2,626263 



C262(62 



I) = D2d'^e2/dX^ , (12a) 



Ads{e3 + 630/2) + ^263 4- ^2(6^ 

= Dsd'^es/dX^ . 



el) + C2e3{ej + el) 
(12b) 
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The term Adsiea + 630/2) in Eq. ( |l2b| ), new with this 
article, results from satisfying the compatibility relations. 
The same equations are obtained, after integrations, on 
using the density of Eqs.(2) plus (3), demanding that F 
be stationary with respect to the displacement, and only 
then using Eqs.(lO). The boundary conditions are 

e2(±oo) = ±\/3e3o/2 , e3(±oo) = -630/2 . (13) 

At A2 — — 2_B|/Ca_ which defines the temperature 
Tbk, the solutions ardj 

62 = (\/363o/2) tanh(«;X) , 63 = -630/2 (14) 

with = 3i?2/(2C2D2); the Ads term and the dilata- 
tional and shear strains all vanish identically. We note 
that Tbk rnay possibly be identified experimentally as 
the temperature where 630 (= 2B2/C2) is three times 
the value at Tc- 

At other temperatures, the strains 61 and eg are not 
zero, and the equations must be solved numerically. We 
define the reduced temperature t — {T — Tq)/{Tc — To); 
then T = 1 at the transition and r = —9 at T = Tbk- 
We also take D3 = D2. 
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32 and 63 as functions of posi- 



tion. The horizontal axis is scaled by k (see Eq. 14) defined 
at Tbk- The parameter Ads is A!2^c — To); the temperatmes 
T ^Tc,T = Tbk and T < Tbk correspond to r = 1, -9 and 
—50 respectively. 

Figure 1 shows 62 and 63 as functions of X at three 
different temperatures, Tc (which is > Tbk), Tbk and 
T < Tbk, as determined from numerical solution of Eqs. 
(12) and (13). The value Ads = A'2{Tc - To) A2(r = 
Tc), which is rather soft, was chosen for display purposes; 
for larger values, 63 remains close to —630/2. One sees 
that 63 + 5630 is > or < for T > Tbk or < Tbk 
respectively. From Eqs. (10) then, the wall region is 
dilated near Tc and compressed below Tbk- 



To estimate the size of these effects, we use data for 
FePd alloys. The data quoted in Ref. | give Tc = 268.6K, 
To — 265K and Tbk = 233K. Combining these with data 
from Rcfs. In], we find that Ads/A2{T = Tc) has a lower 
limit of « 200 based on the smallest observed values for 
A2; we use a more conservative estimate of 400. Between 
Tc and Tbk, the volume change AV/V at the centre of 
the wall reaches a maximum at T « 258 K, with a very 
smah value w lO^^. At T = we find AV/V w -lO^^. 
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FIG. 2. Trajectories in the (63,62) plane for a twin wall 
linking two of the three tetragonal variants (solid circles). The 
parameter Ads and the temperatures are as in Figure 1. 



Figure 2 shows twin-wall trajectories in the (63,62) 
plane. The trajectories bow toward the third variant 
for Tbk < T < Tc, and away forT <Tbk- They shift 
toward the vertical for larger Ads - 

Reference ||, in effect, assumed that ^1 = ylg = 
and so the term Adsies + 630/2) was absent from their 
differential equations. Solution of these equations, done 
only near Tc, gave trajectories which passed close to the 
origin. We point out however that the origin is not the 
cubic state but rather a highly dilated, hig hly sheared 
state with ei — c^q/^/Q, and cq = ± 1/3/8630, as seen 
from Eqs. (10). Reference ||, which also assumed in ef- 
fect that Ai = Aq = 0, proposed a trajectory that is a 
27r/3 circular arc centered at the origin. As evident from 
Figure 2, this is possibly useful only for T ^ Tbk- 

The wall energy W (per unit area) is the energy re- 
quired to form an interface between two variants: 

/oo 
{:Ft^u^n - Th) dX (15) 
-oo 

where Ttwin is the density (11) for the twin- wall solution 
of Eqs. (12) and (13), and 
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= lA2el„ - isae^o + ^^26^0 (16) 

is the density for a single variant. Although not directly 
observable, W has physical content and so we provide 
the following. 
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FIG. 3. Numerical and variational wall energies as func- 
tions of the dimensionless temperature r. The vertical axis 
is scaled by the wall energy at Tbk- From lower to upper, 
the solid curves are results from the solution of the differ- 
ential equations for A^s = (1,10,100) x A'iiTc - To). The 
short-dashed and long-dashed curves are the bounds of Eqs. 
(17) and (18) respectively. 

Figure 3 shows the wall energy W as a function of tem- 
perature, as determined numerically for three different 
values of Ads', it shows also two variational approxima- 
tions which we now obtain. 

If we consider Eq.(14) as a trial fmiction, with k an 
adjustable parameter, we find 

W < eloVsCM^ (17) 

at the optimal k, namely k = 639 -^/SC^/^SZ^ ; the T 
dependence is in 630. From Figure 3, this does very well 
over the temperatures examined, not the least because it 
is an equality at Tbk for all A^s- The coefficients Ads 
and -D3 do not appear because of the form of our trial 
function; explicitly, in obtaining Eq.(17), we did noi drop 
the corresponding terms in the density, and so Eq.(17) is 
valid independent of the magnitude of the dilatational 
and shear energies and strains. Equation (17) seems to 
be an equality also in the limit Ads 00 for|-|all T. 

The wall energy was estimated previouslyQ, using the 
circular trajectory described above. On setting. Ai — 
Aq — Ads = 0, taking — D2, and usingj 63 = 
630 coS(/), 62 — 630 sin0, with ^ < < we find 

W < s/MB^Ih4j27 . (18) 



Unlike Eq.(17), this assumes that Ai = Aq = 0; includ- 
ing dilatation and shear only increases the value relative 
to Eq.(17). Figure 3 shows that Eq.(|T^) is considerably 
poorer than Eq.(17) for the temperatures examined; it 
is a factor of 3 larger at Tc. But the trajectories bend 
toward the Hong-Olson arc at lower T, and so Eq.([l8|) in- 
creases with decreasing T less strongly than Eq.(17) (as 
CgQ^ rather than as Cgg). Indeed, for r < —151, Eq.(p^ is 
actually better than Eq.(17); but this temperature may 
not be accessible since the strain there is ~ 10 times that 
at Tc. For Fe7oPd3o, r = —151 is inaccessible (T < 0). 

In summary, we have found for the first time the solu- 
tion for the C-T twin wall, at all T, in the physical case 
with dilatational and shear energies and strains. The di- 
latation and shear strains, which are localized near the 
interface, change sign at Tbk- The magnitudes (at most 
1% in Fe-Pd alloys) may however be too small to detect. 
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